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We present an ab initio determination of the shear viscosity for the unitary Fermi gas based 
on finite temperature quantum Monte Carlo calculations and the Kubo linear-response formalism. 
The results are confronted with the bound for the shear viscosity originating from hydrodynamic 
fluctuations. Assuming smoothness of the frequency dependent shear viscosity rj[ui), we show that 
the bound is violated in the low temperature regime and the violation occurs simultaneously with 
the onset of the Cooper paring in the system. In order to preserve the hydrodynamic bound in 
QMC r]{Lij) has to possess a sharp structure located in the vicinity of zero frequency which is not 
resolved by an analytic continuation procedure. 



PACS numbers: 03.75. Ss, 05.30.Fk, 05.60.Gg, 51.20.+d 

The unitary Fermi gas (UFG) is expected to be one 
of the strongest correlated systems in nature as it sat- 
urates the unitarity bound for the s-wave cross section 
cr{k) < Aii/k'^ , where k is the relative wave vector of scat- 
tering particles. The strong correlations are responsible 
for a multitude of interesting phenomena, among which 
the existence of the pseudogap regime between the su- 
perfluid state and the normal state [iRsI and the nearly 
ideal hydrodynamic behavior [J-|6[ , are the most surpris- 
ing ones. Therefore the properties of the UFG attract 
enormous attention from several conrmunities including 
atomic physics, nuclear physics, relativistic heavy-ion 
collisions and high-Tc superconductivity, see review pa- 
pers [H3- 

Experimentally the UFG has been realized with 
trapped fermionic atoms by means of Feshbach reso- 
nances 1 1] and currently represents one of the most con- 
trollable quantum systems. On the other hand theoret- 
ical description poses a real challenge due to the lack 
of any small, dimensionless expansion parameter. Fortu- 
nately, for unpolarized UFG the genuine nonperturbative 
quantum Monte Carlo (QMC) techniques are free from 
the infamous fermionic sign problem, and can be used 
to produce reliable predictions, providing a solid bench- 
mark for experimental as well as other theoretical stud- 
ies. Indeed, the method proved its powerful predictive 
capability, providing results which are in very good agree- 
ment with experimental data. In particular, the energy of 
the ground state agrees with experimental results within 
1%[1_2] and the same accuracy was obtained for the equa- 
tion of state for the system being in the normal state [l3| . 
For temperatures around the critical temperature the dis- 
crepancy between experiment and the simulations does 
not exceed 5% [ij]. Other properties well established 
by QMC include the evolution of the pairing gap |15l4l8| 
and the critical temperature for the superfluid-to-normal 



phase transition 19l421|. Last but not least, the QMC 
has been used to determine transport coefhcients of the 



UFG, namely, the shear viscosity [22] and the spin diffu- 
sion cocfhcicnt 23]. 

Determination of transport coefhcients is currently of 
great interest as they provide a verification of various 
hypothetical bounds having their roots directly in quan- 
tum mechanics. The best known example is a conjecture 
formulated by Kovtun, Son, and Starinets (KSS) of the 
existence of a lower bound rj/ s ^ h/iA-nkB) on the ra- 
tio of the shear viscosity 77 to the entropy density s, for 
all fluids [2J]. Due to strong correlations which imply 
short mean free path, the UFG appears as one of the best 
candidates for being a "perfect fluid" (beside the quark- 
gluon plasma created in heavy ion collisions), defined as 
the one with the lowest ratio 77/s. Experimental measure- 
ments for trapped systems yield rj/s ^ 0.5 h/k^ [4-6, 9J], 
while the recent QMC simulations set the minimum value 
to be 77/s w 0.2 h/ kg '22]. Still there are contradicting 
scenarios of how the shear viscosity behaves in the low 
temperature regime. Whereas some authors expect the 
pairing correlations to reduce the viscosity significantly 
and even to decrease it to zero at T = [25, 26], the 
others favor the upturn originating from the phonon con- 
tribution [27*1. 

Recently it was shown by Chafin and Schafer that hy- 
drodynamic fluctuations in a non-relativistic fluid lead to 
the appearance of the minimum in the shear viscosity to 
density ratio rj/n as a function of the temperature [28| . 
In the vicinity of the critical temperature the bound was 
established to be rj/n > 0.3 ft. This result was indepen- 
dently confirmed by Romatschke and Young [29[. While 
the experimental results for trapped systems are consis- 
tent with the hydrodynamic bound, the QMC results for 
uniform systems violate this bound as they predict values 
rj/n < 0.2 ft in the superfluid phase [23] ■ However, to ex- 
tract the shear viscosity within QMC approach one has to 
perform an analytic continuation of the imaginary time 
correlator to real frequencies, which is a highly nontrivial 
procedure. It is known that this step represents an ill- 



posed numerical problem, where statistical fluctuations 
which are tied up with the QMC are greatly enhanced 
during the continuation process. Indeed, it was pointed 
out in the supplement of [2^, and also in [29|, that the 
statistical accuracy of the quantum Monte Carlo signal 
allows to produce the ratio 77/n which is consistent with 
the hydrodynamic bound. 

In this letter we present QMC results for a temperature 
evolution of the shear viscosity of higher accuracy than in 
our previous studies [22| and discuss the errors associated 
with the analytic continuation. Under the assumption of 
the smoothness of frequency dependent shear viscosity 
77(0;), we show that the new results violate the hydro- 
dynamic bound 28|, l29| for 77/n in a low temperature 
regime, while the ratio rj/ s stays above the KSS bound. 
Moreover the hydrodynamic bound is violated for the 
temperatures at which the Cooper pairs are formed in 
the system. We show that the agreement with hydrody- 
namic bound can be achieved only if there exists a low 
frequency sharp peak in the rj{uj) which is overlooked by 
analytic continuation procedure. 

In order to determine the shear viscosity of the UFG we 
employ the QMC technique on the lattice, which provides 
numerically exact results, up to quantifiable systematic 
uncertainties (for details see Ref. 21!]). These simulations 
are very similar to those of Ref. |22|] and therefore we only 
shortly describe the main stages of the computational 
process, focusing mainly on improvements. Henceforth 
we define units: H = m — ks = I- 

We performed simulations using three lattice sizes 
Nx — 8, 10, 12 with corresponding average densities n ^ 
0.08,0.04 and 0.03, respectively. For these lattices the 
systematic errors were estimated to be less than 10—15% 
and they are the most severe in the superfluid phase (for 
detailed discussion in the context of the transport coeffi- 
cients see supplemental materials of Refs. [22', 23]). The 
statistical uncertainties are typically of the order of 1% 
as we collected an ensemble of 10^ uncorrelated samples 
at each temperature. The fact that we have used three 
different lattices allows us to check the stability of the re- 
sults as we approach the the thermodynamic T^ — >■ 00 and 
continuum n — ?► limits, where V = iV^ (the lattice con- 
stant is set to unity). These new results were obtained 
with a new (parallelized) code and for lower densities 
than in Ref. [22], thus closer to the desired limit n -> 0. 
The frequency dependent shear viscosity 77(0;) is obtained 
from the imaginary-time (Euclidean) stress tensor-stress 
tensor correlator 



Gn(g,r) = -(Il(-^)(r)IlLT(0)), 



(1) 



evaluated at zero momentum q = 0, by inversion of the 
relation; 

r f ^ 1 r r ^ coshKr-^/2)] 

Gn(9 = 0,r) = -y^ ,{.)u. 3i,h(../3/2) '^"' ^'^ 



where /3 = 1/T is the inverse of the temperature. Taking 
the limit of zero frequency we obtain the static shear 
viscosity 77 = 7^(0; — > 0). The average in Eq. ([T]) is 
performed over the grand canonical ensemble, 0{t) — 
^r(H-t.N)Q^-T{H-t,N)^ fj jg ^j^g Hamiltonian of the sys- 
tem, /i is the chemical potential, and N is the particle 
number operator. For evaluation of the shear viscosity 
it is sufficient to use the kinetic part of the stress tensor 
only (see e.g. Ref. [30J): 
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The correlator Gu{q = 0, r) was evaluated for at least 61 
points in imaginary time t, uniformly distributed in the 
interval [0, /?]. Since the correlator is symmetric with re- 
spect to /?/2 point, we restrict the inversion procedure to 
the interval t S [0,/?/2], and apply the symmetrization: 
Gn{T) ^— [Gn{T) + Gn(/3 — t)]/2 in order to decrease 
statistical fluctuations. 

The inversion of the Eq. ([2|) represents numerically an 
ill-posed problem. However, by using known theoretical 
constraints the solution of this problem becomes easier 
to determine. Besides the non-negativity of the shear 
viscosity ri{uj) ^ 0, the sum rule 



dio 



riiuj) 



C 



Ibn^/uJ 



(4) 



and the asymptotic tail behavior 77 (w — > 00) = 
C/{157r^/uJ) have been used as a priori information. The 
Tan contact density C and the energy density e are 
obtained consistently within the same simulation, and 
are in agreement with results reported in previous stud- 
ies 21J, |23|. To perform the inversion we applied the 
methodology which combines two complementary meth- 
ods: singular value decomposition (SVD) and maximum 
entropy method (MEM), both described in Ref. [3l|. The 
inversion procedure consists of two steps: i) the SVD 
method provides us with the projection of the solution 
onto the subspace, where the inverse problem is well- 
posed 77, ii) the inversion is performed using the self- 
consistent MEM and we impose the SVD solution as ad- 
ditional external constraint i.e. ^[77] = 77, where P stands 
for the projection operator onto SVD subspace. In or- 
der to estimate the stability of the method with respect 
to the algorithm parameters, the "bootstrap" strategy 
was applied. For more details see supplemental materi- 
als of [gi, Q . 

An important ingredient of the self-consistent MEM 
is an appropriately chosen class of a priori models for 
the solution. We assume that the frequency dependent 
shear viscosity 77(0;) has Lorentzian-like structure at low 
frequencies, smoothly evolving into the asymptotic tail 



behavior: 



C 



where 



and 
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The parameters {/i, 7, 0,02,02} describe adraissible de- 
grees of freedom of the model and are adjusted automat- 
ically in the self-consistent manner. This choice of model 
space is reasonable as it is compatible with the low fre- 
quency behavior of the 77(0;) predicted by hydrodynamic 
fluctuations formalism 28, 2!J] and with results obtained 
within Luttinger-Ward theory |30j. Finally, we assume 
that there is no sharp structure in the shear viscosity 
77(0;) in the low frequency limit. This assumption is of 
crucial importance, as it is well known, that sharp struc- 
tures can be overlooked by numerical analytic continu- 
ation [3l[. Clearly, note that adding e(5(w) to the shear 
viscosity r\(y>) corresponds to increasing the imaginary 
time correlator by a constant value AGn — e/7r/3. 

In Fig.[l] the dimensionless frequency dependent shear 
viscosity 'q(uj)/n is shown for two selected temperatures 
obtained for 10^ lattice: close to the critical temperature 
of supcrfluid to normal phase transition Tc — Q.\b{\) £p 
and in the normal phase T = 0.30 e^, where e p — ^kp — 
i(37r^n)^/"^ is the Fermi energy of the non-interacting gas. 
In the normal phase T > 0.25 e^, one may observe that 
r){ijj) possesses a Lorentzian-like structure at low frequen- 
cies with a maximum in the vicinity of zero frequency. 
For temperatures below T < 0.25 e^. the suppression of 
77(0;) develops for the zero frequency. We associate this 
depletion with the onset of pairing correlations, which 
has been identified to exist in UFG above the critical 
temperature up to temperatures T* = 0.20 — 0.25 £p., 
commonly referred to as pseudogap regime 13, [l^ [23[ . 

In Fig. [21 the dimensionless static shear viscosity rj/n is 
shown as a function of dimensionless temperature T/sp. 
Results obtained for different lattices exhibit reasonable 
agreement, even though they are affected by significant 
uncertainty of the order of 10% — 20%. The tempera- 
ture dependence of the shear viscosity can be resolved, 
which is clearly visible when a smooth approximation 
of QMC results is performed, while the amount of scat- 
tering of the QMC points can be used to estimate fi- 
nal uncertainty of the viscosity computation - yellow 
band in Fig. [2l Moreover, the new results are consis- 
tent with previous estimate of the viscosity coefficient of 
Ref. J22| . In the low temperature regime no upturn orig- 
inating from phonon contribution is observed, however 




FIG. 1: (Color online) The frequency dependent shear vis- 
cosity obtained from the QMC calculations for 10^ lattice for 
two different temperatures: upper panel T = 0.16 ep (close to 
the critical temperature) and lower panel T — 0.30 ep (normal 
phase). The gray bands show uncertainty of the reconstruc- 
tion. In the insets, the corresponding QMC correlators (blue 
points) are presented. The lines show quality of the repro- 
duction of the QMC data. 



it does not imply that there is no upturn at all. Pre- 
cisely, the boxes used in the simulations are not capa- 
ble properly capturing the phonon related physics. It 
becomes clear if we note that the phonon occupation 
probability nph(fc) ex fc^/[exp(cfc/T) — 1] (where c is 
the sound velocity) has a peak at momenta of the or- 
der of fc « 2.2 kpT/EF, which for T/sp = 0.10 •••0.15 
is smaller than the smallest lattice momentum 2tt/Nx. 
Thus infrared momentum cutoff 2t:/Nx removes almost 
entirely the pure phonon contribution. To take into ac- 
count this source of the systematic error we modified 
the estimated error of the shear viscosity computation 
in such a way to encapsulate the pure phonon contribu- 
tion of Ref. 



27[. Phonons appear to contribute signifi- 



cantly only well below the critical temperature Tc. Above 
the critical temperature, when the system is in a normal 
state and the superfluid order parameter is vanishing the 
calculation approach adopted in Ref. '27'] is not valid 
anymore. This approach is based on an effective theory 



1.2 
1 
0.8 
0.6 
0.4 
0.2 






* * 



QMC average 

classical hydro. +therm.fluct. 

phonons contribution 

QMC (2012) 

N^=8, n=0.06 

N^=10, n=0.04 

N.=12, n=0.03 



0.1 



0.2 T 



0.3 



0.4 



0.5 



T/fip 



FIG. 2; (Color online) The dimensionless static shear vis- 
cosity ri/n as a function of temperature for an 8^ lattice - 
solid (red) circles, 10^ lattice - (blue) squares and 12'^ lat- 
tice - (green) diamonds. Vertical black dotted lines indicate 
the critical temperature of superfluid to normal phase transi- 
tion for Tc = 0.15 Ep and the onset of Cooper pair formation: 
T* « Q.22 Ep, respectively. The black line indicates smooth 
approximation of the QMC results and the (yellow) band re- 
flects uncertainty of the viscosity computation. The (purple) 
asterisk shows the shear viscosity from the QMC calculations 
of Ref. ^21 ■ The phonon contribution to the viscosity, eval- 
uated according to Ref. ^3 is shown as a dashed (brown) 
line and it is accounted for in the (yellow) band. The hydro- 
dynamic bound of Ref. [2y] in plotted by black dashed line. 
The inset shows ratio of the shear viscosity to entropy den- 
sity {rj/s) as a function of T/Ep, obtained from the smooth 
approximation of the QMC results. The black dashed line 
indicates KSS bound: l/47r. 



for the phonon field alone, valid at very low temperatures 
only, where Cooper pair breaking and pseudogap physics 
can be ignored. 

Comparing our ah initio results with the hydrodynamic 
bound we find that as long as the system is in the normal 
phase (for T > T*) the static shear viscosity is located 
above the bound. The results obtained for all three lat- 
tices indicate, that the hydrodynamic bound is violated 
when the system enters the pseudogap regime, which oc- 
curs above the critical temperature, at T* w 0.25 e p. In- 
terestingly, the corresponding ratio 77/s stays above the 
KSS bound rj/s ^ l/47r « 0.08 for all considered temper- 
atures, with the minimal value (?7/s)jjjijj « 0.2, located 
in the vicinity of the critical temperature. For the tem- 
perature evolution of the entropy we have used the MIT 
results of Ref. [331 j ^^ o^^ results reproduce experimental 
data with reasonable accuracy, see |l4l22|. 

The presented methodology, which produces the most 
probable results, together with the smoothness assump- 
tion of r]{uj) clearly allows for the hydrodynamic bound 
to be violated in the low temperature regime for T <T* . 
The hydrodynamics bound was obtained by adding the 
contribution due to thermal fluctuations to an undefined 



classical value of the viscosity 28|, |29[ and with an am- 
plitude depending on a cutoff parameter, which can be 
estimated only by order of magnitude, and thus a fac- 
tor of order unity cannot be refuted. These estimates 
of a possible lower bound for viscosity are purely classi- 
cal in character and they show that the role of classical 
thermal fluctuations could become important if the bare 
viscosity, evaluated in the absence of such thermal fluc- 
tuations, is small. The bare viscosity can be evaluated 
in a classical kinetic approach, but to what extent such 
an estimate is valid in the case of a unitary Fermi gas 
remains an open question. A unitary Fermi gas is not 
dilute in the sense that only two-body collisions are im- 
portant, since am?^^ — 0(1), where a ~ 4:7r/kjp is the in- 
teraction cross section between fermions at the Fermi sur- 
face and n — kp/Sir'^. Even though the role of very low 
momentum phonons at very low temperature cannot be 
properly included in QMC within the presently tractable 
lattice sizes, still all other effects (pairing, pseudogap, 
fermion collisions) are accounted for. As seen from Fig. [5] 
phonons show a significant contribution only for temper- 
atures below the critical temperature T^. It implies that 
even if the phonons were included, the viscosity would 
have remained below the classical hydrodynamic bound 
23,123 up to temperatures smaller than T « O.le^. 

Even though we have no reasons to expect that any 
alleged unaccounted for thermal fluctuations of the kind 
advocated in Refs. [28|, |29| might appear in the pseudo- 
gap regime below the temperature T* , we will force such 
a contribution into the spectral function to investigate 
whether our QMC data are consisted with this classical 
hydrodynamical bound. We have performed three tests 
using data obtained for 10'^ lattice at the critical tem- 
perature. Hydrodynamic fluctuations '28|, |29|, suggest 
non-analytic contributions for small frequencies of the 
form: 



T^iuj) = r;(0) - V^T 
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(8) 



where the momentum diflusion constant D,^ is positive. 
Thus, 77(1^) is predicted to be monotonically decreasing 
function for small frequencies. In the first test we forced 
the MEM solution to be consistent with this assumption 
by replacing the Lorentzian-like structure at low frequen- 
cies (O by the expression: 



Ciuj,{A,B})^A-By^, 



(9) 



where A> 0, B > 0. This modification implies that the 
space for default models contains only solutions which are 
monotonically decreasing and explicitly consistent with 
predictions taking into account hydrodynamic fluctua- 
tions. In addition we have released the constraint pro- 
vided by SVD. The resulting frequency dependent vis- 
cosity is shown with a (red) solid line in Fig. [3] In the 
second test we fitted the modified model to the imaginary 
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FIG. 3: (Color online) The frequency dependent shear vis- 
cosity for the temperature T = 0.16 £^ obtained by three 
different methodologies of the reconstruction as described in 
the text. The (purple) asterisk indicates the hydrodynamic 
bound obtained in Refs. [2S.i29j|. 



time correlator, in a similar way as it has been done in 
Ref. [29[ . This method typically violates the sum rule ^ 
and the resulting r]{ui) is shown with (blue) dot-dashed 
line in Fig. |3l Finally we fitted the QMC correlator to 
the model proposed by Enss et al. in the paper [30| 



ry(w) 



Wt„ 



C WT^(1 + LUTjj) 



1 + (u;t,,)2 IStt^^ 1 + {ujTr,y 



(10) 



where the only free parameter is the viscous transport 
scattering time r^, as the total Drude weight W is fixed 
by the sum rule dH), see (red) dotted line in Fig. [31 All 
three approaches are able to reproduce the imaginary 
time correlator within its error bars. These simple tests 
result in a static shear viscosity below the hydrodynamic 
bound. In the case of MEM with the model ([9]) we find 
that the i]/n typically is located below the bound. How- 
ever we find settings of the MEM algorithm which pro- 
duce results above the bound as refiected in error band 
depicted in the Fig. [3l Similar results we find for other 
temperatures. 

Based on these tests we conclude that the QMC re- 
sults are compatible with the hydrodynamic bound only 
if we assume the existence of a very narrow peak located 
around a; w 0. The classical hydrodynamic fiuctuation 
formalism predicts this type of peak, which has its ori- 
gin in non-analytic contributions to the shear viscosity. 
In Rcf. [28] it was estimated that the width of the non- 
analytic structure can be of the order of Aw ~ 0.3r(n/77) 
and at the critical temperature it gives Auj/sp ~ 0.15, 
which is in agreement with our first test. However, we 
have to point out several important aspects of our results. 
i) We do not observe any abrupt change in the frequency 
dependence shear viscosity profile as we decrease the tem- 
perature, which is a typical signal of losing the resolution. 



Instead, the profile we construct evolves smoothly be- 
tween two cases presented in the Fig. [1] ii) Our QMC 
simulations admit the existence of the pseudogap regime 
for the UFG 
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23l |. Other theoretical approaches 
which admit the presence of the pseudogap predict the 
shear viscosity to be significantly suppressed as in the sys- 
tem a fermionic gap opens at T* and they clearly violate 
the hydrodynamic bound |25l . |26| . Similar behavior was 
predicted for the spin diffusion coefficient [33]. Our ab 
initio results are qualitatively consistent with these ob- 
servations and support scenarios provided by these the- 
ories. Hi) We find it difficult to rationalize the existence 
of a non-analytical classical hydrodynamical correction 
to the spectral function deep in the quantum regime. 

In summary, we have presented the temperature evolu- 
tion of the shear viscosity obtained within fully ab initio 
simulations. Assuming the smoothness of the ri{uj) the 
results violate the hydrodynamic bound for rj/n in a low 
temperature regime, while the ratio r]/s stays above the 
KSS bound for all studied temperatures. Violation of the 
hydrodynamic bound occurs simultaneously with occur- 
rence of the Cooper pairs in the system. The accuracy of 
the QMC signal allows a result with the static viscosity 
being above the hydrodynamic bound only if we assume 
an existence of a sharp structure in the ?y(a;) located in 
vicinity of zero frequency. 
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